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Introduction 

Considerable progess has been made in recent years in chemical functionalization of fullerene 
molecules [see, for example: Taylor 93, and Diederich 96]. In some cases, the predominant 
reaction products are different from those obtained (using the same reactants) from polycyclic 
aromatic hydrocarbons (PAHs). One such example is the cycloaddition of o-benzyne to C60. It 
is well established that benzyne (Fig. 1) adds across one of the rings in naphthalene, 
anthracene and other PAHs [Hoffmann 67] forming the [2+4] cycloaddition product 
(benzobicyclo[2.2.2.]-octatriene with naphthalene and triptycene with anthracene (Fig. 2)). 
However, Hoke et al [Hoke 92] 



Fig. 1 Space-filling model of benzyne. Gray atoms are carbon and white atoms are hydrogen. 
Fig. 2 Triptycene, the product of benzyne addition to anthracene. 

demonstrated that the only reaction path for o-benzyne with C 69 leads to the [2+2] 
cycloaddition product in which benzyne adds across one of the interpentagonal bonds 
(forming a cyclobutene ring in the process). Either reaction product results in a loss of 
aromaticity and distortion of the PAH or fullerene substrate, and in a loss of strain in the 
benzyne It is not clear, however, why different products are preferred in these cases. Fig. 3-5 
are the [2+2] and [2+4] products for addition of benzyne to naphthalene and C60, respectively. 
There are no experimental data available for the addition of benzyne to large fullerenes or 
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carbon nanotubes. However, if multiple benzyne molecules could be added to nanotubes at 
controlled sites, the resulting products could resemble molecular paddle wheels or gears and 
might be useful components for nanotechnology machines (this concept is investigated in a 
companion paper [Han 96]). 



Fig. 3 Products of the benzyne-naphthalene cycloaddition reaction, (a) [2+2] product, (b) [2+4] 
product. 



Fig. 4 [2+2] product of the benzyne-C60 cycloaddition reaction (two views). 
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Fig. 5 [2+4] product of the benzyne-C60 cycloaddition reaction (two views). 

In the current paper, we consider the stability of benzyne-nanotube adducts and the ability of 
Brenner’s [Brenner 90] potential energy model to describe the structure and stability of these 
adducts. The Brenner potential has been widely used for describing diamondoid and graphitic 
carbon. Recently it has also been used for molecular mechanics and molecular dynamics 
simulations of fullerenes and nanotubes. However, it has not been tested for the case of 
functionalized fullerenes (especially with highly strained geometries). We use the Brenner 
potential for our companion nanogear simulations [Han 96] and believe that it should be 
calibrated to insure that those simulations are physically reasonable. In the present work, 
Density Functional theory (DFT) calculations are used to determine the preferred geometric 
structures and energetics for this calibration. The DFT method is a kind of ab initio quantum 
chemistry method for determining the electronic structure of molecules. For a given basis set 
expansion, it is comparable in accuracy to the MP2 method (better than Hartree Fock, but less 
accurate than more extensive electron correlation methods such as MP4 or CCSD). However 
for systems with large numbers of basis functions it more efficient than any other methods that 
include electron correlation effects. 

In this presentation we show the results of DFT calculations for the reaction of benzyne with 
naphthalene, C60, and nanotube models. We compare energies for [2+2] and [2+4] 
cycloaddition products. The preferred products for the naphthalene and C60 reactions have 
been determined by experiment and, thus, these cases serve as a validation of our quantum 
chemical approach. We also compare the DFT and Brenner potential results. Finally we can 
predict the likelihood of reaction between benzyne and nanotubes. 

Methods 

We have carried out calculations for the addition of benzyne to naphthalene and C60 where 
the preferred products have been established by experiment. Both the [2+2] and [2+4] 
cycloaddition products have been considered in these cases. Then we describe calculations 
for both products of the benzyne-nanotube reaction. The calculations were carried out on IBM 
RS6000 model 590 workstations using a MULLIKEN_2.0 [Rice 96]. The STO-3G basis set 
was used for the geometry determinations, unless otherwise specified. The DFT calculations 
were carried out using a hybrid functional containing Hartree Fock exchange and electron 
correlation terms. Specifically, we use the 3-parameter empirical formulation devised by Becke 
[Becke 93] that is usually denoted as B3LYP. The parameters were determined by Becke to 
best reproduce atomization energies and ionization potentials for a set of 56 small molecules 
and is widely regarded as the best all-purpose DFT functional. For some of the cases studied, 
we carried out calculations using the larger 4-31G basis set at the STO-3G/DFT geometries. In 
addition, for the naphthalene-benzyne system we also obtained geometries and energies 
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using the 6-31 G* basis set. These basis sets are standard sets and widely used. STO-3G is a 
minimal basis set having two s atomic orbitals and one set of p orbitals for each carbon atom 
and a single s orbital for hydrogen. 4-31 G is a near double zeta basis set containing three s 
and two sets of p orbitals for carbon and two s orbitals for hydrogen. 6-31 G* is similar to 
4-31 G, but it also includes one set of d orbitals for each carbon. Each orbital is comprised of a 
linear combination of up to 6 gaussian functions. For each carbon atom these basis sets result 
in 5 (STO-3G), 9 (4-31 G), and 14 (6-31 G*) atomic orbitals. 

For the nanotube calculations, a segment of a single-walled [14,0] nanotube was generated 
(see [Han 96]). This tube has 14 benzene rings around the circumfrence and a diameter of 
10.9 A. The segment used in the present work is only 3 rings long, but this is adequate to 
describe the distortions that occur at the site where the benzyne is joined. In order to make the 
quantum chemistry calculations tractable, a further simplification had to be made: we used 
only a semicircular band to represent the nanotube, but constrained the carbon atoms at the 
two ends of the strip to remain 10.9 A apart. This nanotube model with the benzyne adduct is 
given in Fig. 6 and 7. The benzyne carbon atoms are shown in red for easier visibility. The 
blue atoms are the end carbons subject to the distance constraint. Keeping the two ends of 
the strip 10.9 A apart forces it to maintain the appropriate curvature. This nanotube model 
(without the benzyne adduct) has the chemical formula C64H22. 



Fig. 6 Model used for benzyne-nanotube reaction: [2+2] product. Red atoms are benzyne 
carbons; blue atoms are nanotube carbons subject to distance constraints. 
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Fig. 7 Model used for benzyne-nanotube reaction: [2+4] product. 

Results 

The optimized structures for the benzyne-naphthalene adducts are shown in Fig. 3a for the 
[2+2] cycloaddition product and in Fig. 3b for the [2+4] cycloaddition product. The benzyne is 
brought down from above and causes the naphthalene to fold. In the former case this folding 
is along the interior C-C bond and the benzyne attaches by a 4-member ring (cyclobutene). In 
the latter case one of the naphthalene rings folds and the attachment is achieved by two 
6-member cyclohexadiene rings. Both products have C2v symmetry. The relative energies of 
the two product molecules are given in Table 1 . It is seen that the [2+4] product is 
considerably more strongly bound than the [2+2] product, consistent with the experimental 
observations. It is also seen that the binding energies are only moderately sensitive to the 
basis set used, and the energy difference between the products is essentially independent of 
basis set. 


Table 1. Energy of Reaction (in kcal/mol) for Benzyne 

Cycloaddition 


j Reactant 

Product [ 

STO-3G 

4-31 G 

6-31 G* 

Brenner 

| Naphthalene 

[2+4] | 

-67 

-65 

-66 

-96 


[2+2] 

-4 

-8 

-12 

+5 

| C60 

[2+4] 

-25 

+10 


-107 


[2+2] 

-108 

[ -68 


-109 

jNanotube[14,0] 

[2+4] [ 

-88 



-108 


| [2+2] :[ 

-56 



-83 


For the [2+2] product, the interior naphthalene bond length increases from 1 .41 to 1 .60 A and 
the benzyne "triple" bond increases in length from 1 .22 to 1 .38 A. The bonds connecting 
benzyne to naphthalene are 1 .53 A long. The 4-membered ring is a planar trapezoid with 
angles of 95 and 85°. The dihedral angle describing the folding of the naphthalene changes by 
49.4° between reactant and product (from 180 to 130.6°). Mulliken population analysis shows 
a modest increase of negative charge at the naphthalene connection points from -.004 to 
-.127 and a change in charge at the benzyne connection points from -.025 to +.063. For the 
[2+4] product, the connecting bonds are 1 .56 A long. The benzyne sits 1 .44 A above the 
original naphthalene ring plane. However the resulting adduct has 3 nearly equivalent 
cyclohexadiene rings and the original ring of the naphthalene is bent by 60° from planarity. 

These naphthalene-benzyne adducts were subjected to various distortions to determine the 
stiffness of these molecules as an indicator of how well they might function as teeth in 
nanogears [see Han 96]. Tilting and twisting distortions refer to the position of the benzyne 
relative to the underlying naphthalene molecule. Tilting is the distortion largely responsible for 
slip in the nanogears [Han 96]. For a 10° tilt, the energy of the [2+2] product increased by 1 .32 
kcal/mol (5.56 kcal/mol for a 20° tilt) and the energy of the [2+4] product increased by 3.94 
kcal/mol. Thus, it is seen that the [2+4] product is somewhat stiffer than the [2+2] product. A 
10° twist increased the energy of the [2+2] product by 0.82 kcal/mol. We also pulled the 
benzyne ring away from the naphthalene. For the [2+2] product, when the benzyne is pulled 
off by 0.4 A the resulting energy is 14.5 kcal higher than the asymptote, indicating the 
existence of an energy barrier. At this geometry, the naphthalene had almost returned to 
planarity (as can be seen in Fig. 8). 
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Fig. 8 Partially dissociated benzyne-naphthalene [2+2] adduct. Benzyne pulled back 0.4 A 
from equilibrium value. 

The geometries of the [2+2] and [2+4] benzyne-C60 adducts are shown in Fig. 4 and 5, 
respectively. Two view are given for each product. The relative binding energies are given in 
Table 1. In this case, the [2+2] product is considerably more strongly bound than the [2+4] 
product, also consistent with the experimental observations. 

Finally, the benzyne-nanotube adducts are shown in Fig. 6 and 7 and their relative energies 
given in Table 1 . In contrast to the other systems studied, both products are strongly bound. 
However, the [2+4] product is more stable than the [2+2]. 

Examination of Table 1 shows that the Brenner potential qualitatively describes these products 
fairly well. The main failing is that it predicts the [2+4] C60-benzyne product to be almost as 
stable as the [2+2]. It also predicts the naphthalene and nanotube products to be 20-30% too 
strongly bound compared to the DFT results. However, that might be due to limitations in the 
quantum chemistry calculations. Not shown in the table, however, is the fact that the Brenner 
potential predicts the existence of other stable C60-benzyne adducts that are not observed in 
experiment. For example, a [2+2] product with the benzyne attached to a C-C bond joining 5- 
and 6-membered rings is slightly more stable than the one joining two 6-membered rings (the 
one observed in experiment and quantum chemistry studies). We are currently carrying out 
quantum chemistry calculations to look for this alternate product. 

Conclusions 

We have carried out quantum chemistry calculations using the DFT method to study possible 
products of the reaction of benzyne with PAHs, fullerenes and nanotubes. The results for 
naphthalene and C60 are in good agreement with experimental observations and confirm that 
the benzyne-PAH reaction leads to a [2+4] cycloaddition product and the benzyne-fullerene 
reaction leads to a [2+2] cycloaddition product. These two products for the benzyner- 
nanotube reaction are strongly bound and should both be stable. The [2+4] product is more 
stable than the [2+2] product. Qualitatively the Brenner potential agrees with the quantum 
chemistry calculations regarding these general conclusions. However, it appears to not 
discriminate between the products of the benzyne-C60 reaction and thus does not predict the 
product selectivity seen in the experiments. 
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